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Convective counterparts of variants of the nonlinear Fisher equation which describes reaction 
diffusion systems in population dynamics are studied with the help of an analytic prescription and 
shown to lead to interesting consequences for the evolution of population densities. The initial value 
■ problem is solved explicitly for some cases and for others it is shown how to find traveling wave 

solutions analytically. The effect of adding diffusion to the convective equations is first studied 
through exact analysis through a piecewise linear representation of the nonlinearity. Using an 
appropriate small parameter suggested by that analysis, a perturbative treatment is developed to 
» ! ' treat the case in which the convective evolution is augmented by a small amount of diffusion. 
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I. INTRODUCTION 



The Fisher equation, proposed originally to describe the spread of an advantageous gene in a population 
has found a great deal of use in mathematical ecology, directly 0, 0, 0, as well as indirectly, i.e., in modified 
M-J | forms. Examples of modifications are the incorporation of internal states such as those signifyin g th e presence 
Oh . or absence of infection 0, 111, the introduction of temporal nonlocality in the diffusive terms |gL Il0|. of spatial 
nonloc ality in the competition interaction terms |lllll2j |. and the addition of convective terms signifying 'wind effects' 
Hia, Uj> UH HH 03 Most of the modifications lead to equations which, as in the case of the original Fisher equation, 
are not soluble analytically and require approximate or numerical methods for their analysis. There is, however, one 
modification |l5l |. that leads to an analytic solution via a trivial transformation, and to the possibility of the extraction 
of interesting information of potential use to topics such as bacterial population dynamics. That modification consists 
of the replacement of the diffusive term in the Fisher equation 

du ( u\ r^d 2 u 

— = au 1 )+D—- (1) 

dt V KJ dx 2 y ' 

o ■ 

ly-j , by a convective term. As is well-known, in the original Fisher equation, u represents a concentration or density, a the 
growth rate, K the so called 'carrying capacity', and D the diffusion constant. In spite of the fact that purely convective 
, (non-diffusive) nonlinear partial differential equations are much less sophisticated than their diffusive counterparts, 
there are at least two reasons why their study is important. There are clear physical situations [la. Il4. flq in which 
q \ wind effects which add the convective term to a diffusive equation are present. These arise, for instance, when one 
changes one's reference frame to that of a moving mask in bacterial population dynamics experiments |l3| . Because 
the convective term comes from the macroscopic motion of the mask under experimental control whereas the diffusive 
term comes from the microscopic motion of the bacteria, it is possible to arrange the system experimentally so that 
the diffusion effects are relatively unimportant. This may be achieved for instance, by using viscous environments, or 
by genetically engineering the (bacterial) population. In such a case one arrives at the analysis of a largely convective 
(neglig ibly diffusive) nonlinear evolution which may be investigated in zeroth order as a purely convective equation 
|l5lll7| |. This fact, that experimental realization of the convective nonlinear equation is indeed possible, provides one 
motivation for these studies. The other, as explained in Ref. 15], is that it is possible to develop perturbation schemes 
starting from the soluble convective equation to incorporate diffusive effects. In Sections ITU and ITTT1 we develop the 
theory for arbitrary nonlinearities in the absence of the diffusion constant, and in Section llVI we focus on the effect 
of adding the diffusive term. 

Recently, a simple prescription for generalizing the analysis of Ref. was given by one of the present authors 
|17| to arbitrary nonlinearities in the convective equation. The equation considered there is 

du du „. . 

_+„_ = «*/(„), (2) 

where v is the medium or 'wind' speed, a is a growth rate and f(u) is a dimensionless arbitrary nonlinearity. The 
prescription provided in that analysis gives the full time and space evolution of u(x, t) in terms of its initial distribution 
u(x,0) = uo(x), through the relation 

u(x,t) = G- 1 [G[u (x-vt)]-at] (3) 
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where the function G{u) is obtained from the nonlinearity /(it) by the integration [l^ 



We have used a slightly changed notation here relative to ref. 0] . Our functions G(u) and f(u) are both dimensionless 
in contrast to the dimensioned counterparts g(u) and F(u) in ref. 

In this paper we present several applications of the prescription and report a number of additional results from 
the nonlinear convective equation (0). In Section [HJ we examine features of the solutions obtained exactly from the 
prescription for several cases of the nonlinearity including an interesting 'pyramid effect' that occurs for nonlinearity 
functions f(u) with multiple zeros. In Section lTTTl we develop a general analysis for finding traveling wave solutions, and 
to study characteristics of the tails and shoulders of such solutions. Applications of this analysis for the nonlinearities 
studied in Section ^ are also made in Section lllll In Section II VI we analyze in detail the perturbative problem 
obtained by adding a small amount of diffusion to the nonlinear convective equation. Conclusions are presented in 
Section |V| 



II. APPLICATIONS OF THE PRESCRIPTION FOR THE SOLUTION OF THE INITIAL VALUE 

PROBLEM 

Applications of the prescription @ for the solution of the arbitrary initial value problem were indicated in ref. [l7| 
for three cases: the ordinary logistic nonlinearity for which f(u) has a single zero at u = K, the Nagumo nonlinearity 
for which f(u) has two zeros, and the trigonometric nonlinearity for which f(u) has multiple zeros. We will display 
new results for the latter two, and analyze other generalizations of the logistic nonlinearity. In some of these cases 
(subsections A,D and E below) u(x,t) is presented in an implicit form, whereas in the others (subsections B and C) 
the form is given explicitly. 



A. Allee effect from Nagumo nonlinearity 

In J2J let f{u) be given by 

«•)-(> "*)(*-*) < 5 > 

where A < K. Such an f(u) represents the so called 'Allee effect'. This effect 0] is associated with the existence 
of an additional zero (fixed point) in the nonlinearity relative to the logistic case. Unlike in the logistic case, the 
zero-it solution is stable here. If u is small initially, it is attracted to the vanishing value; if large, it is attracted to 
the non-zero value u = K . The demarcation point is the additional zero introduced in this case, viz. u = A. The 
physical origin of the Allee effect in population dynamics is the possible increase of survival fitness as a function of 
population size for low values of the latter. Existence of other members of the species may induce individuals to live 
longer whereas low densities may, through loneliness, lead to extinction. There is ample evidence for such an effeect 
in nature [l^, |2fJ, |21| .To construct G(u) from (@J is straightforward, 



_L f {K-u) A u K - ' 

A{K^~A) n ^ V^~M 



~ TTrT — IT m ( i~ a\k ). W 



but to invert G{u) is not. In Fig. Q the nonlinearity is displayed along with the time evolution for a certain initial 
condition. The Allee effect causes a decay to zero of parts of the solution, while other parts evolve to the saturation 
value K. Thus, an initial condition with u values belonging to both regions, will show a growth for those points in 
space such that uq(x) > A, while those points such that uq(x) < A will decay to zero. The initial condition in this 
case evolves into a step profile of width given by the separation in space of the two points where uq(x) = A. 



B. Trigonometric nonlinearity and the Pyramid Effect 

What would happen if the nonlinearity function uf(u) had multiple zeros rather than three as in the Nagumo 
case? The physical origin of such behavior may lie in nonmonotonic response of survival fitness to population size. 
Complex biological and sociological interactions among members of the species could exist on the basis of the coupling 
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FIG. 1: The Nagumo nonlinearity of Eq. (2J and its effects on evolution. The nonlinearity is plotted in (a). Shown in (b) is 
an example of the evolution from an initial condition sech(x/a) (shown by the dotted line) whose profile belongs both to the 
region where < u(x) < A and where A < u(x) < K. The portion in the first (lower) region decays while that in the second 
(upper) region grows to the saturation value u(x) — K. Here, a is an arbitrary length. 



of mating instincts with mere size effects. Such interactions could, in principle, make the survival fitness optimum 
at several, rather than a single, values of the population size. In order to address such situations, extensions of the 
normal Alice considerations are necessary. For this purpose, the discussion in A above may be generalized as follows. 
For a given nonlinearity with an arbitrary number of zeros, uf(u) can be sectioned among each subsequent unstable 
zero. The dynamics of each point of the initial profile will evolves according to which of these sections it belongs 
initially. In the Nagumo case there are only two such regions: the first region is between u = and u = A, while the 
second region is for u > A. All the points of the initial condition such that uq(x) < A decay to zero, while all those 
points belonging to the second region grow to the saturation value u = K . But when there is more than one unstable 
zeros, as in the case of a trigonometric nonlinearity, the propagating solution may evolve into a pyramidal profile. To 
show this we choose a sinusoidal nonlinearity given by 



/(«) 



K sin (at) 



(7) 



where we took the opportunity to change the stability of u — compared with the Nagumo case presented in A. From 
G(u) is given by 



and the exact solution at all time is 
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(9) 



Let us consider an initial condition such that < uo(x) < 2mK where m is a positive integer. All those initial 
points of the profile that belong to the region delimited by (2m — 2) < u(x) < 2mK can be described by the 
following dynamics. The solution grows from u = (2m — 2)K to the value u — (2m — l)K and decays to it from 
u = 2mK following the evolution of an initial condition with compact support in that region. As extensively studied 
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in the case of the logistic nonlinearity |15(, any initial condition with a bounded domain evolves into a step profile 
of width given by its initial value. In the sinusoidal nonlinearity, from each region a step profile with width given 
by [2mK/ir — u (x)] — 9 [(2m — 2)K/tt — uq(x)] eventually emerges. Figure (0 shows an example for an initial 
condition that overlays three such regions. It is clear from the figure that in the bottom region the profile behaves as 
if the other regions did not exist. 




FIG. 2: The 'pyramid effect' for a nonlinearity with multiple zeros. The sinusoidal nonlinearity of Eq. (0 is plotted in (a). The 
horizontal lines at u/K — 1,3,5 depict levels of separation of the regions between stable zones. The corresponding evolution 
is plotted in (b) for the initial profile uo(x) — 5.9 K sech(a;/cr), where a is an arbitrary length. The initial profile is shown as a 
dotted line. The profile moved rightwards as shown. 



C. Multi-individual struggle for environment resources: linear growth with power saturation 



One can envisage physical situations in which one encounters a simple generalization of the logistic nonlinearity 
that retains the linear behavior near u = but accentuates the nonlinearity for larger values of u. This would arise, 
e.g., if the elemental struggle of the individuals for resources involves not two but multiple individuals. The struggle 
would then not be binary in nature. To address this possibility, we consider 



f(u) = [1 - (u/KY 



(10) 



The saturation effect already seen in the logistic case, becomes more abrupt as n increases. Pictorial representations 
for several cases of n are provided in Fig. © . Application of the prescription Q produces the explicit form for G(u) 



G(u) = - In 



1 - (u/K) 7 



(u/K) 1 



(11) 



Its inverse can be written analytically. An arbitrary initial condition uq(x) = u(x,0) therefore via prescription @ 
evolves according to 

i*o (x — vt) 



u(x, t) 



*) + e- 



1/r. 



(12) 



D. Nonlinear growth with power saturation as in bisexual reproduction 



Interesting consequences arise if the nonlinearity cannot be approximated by a linear term near u — 0, such as when 
the growth of the species is nonlinear. A physically relevant example is the representation of bisexual reproduction 
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FIG. 3: The nonlinearity with power saturation and its effect on the time evolution. The nonlinearity with f(u) given by HUH 
is shown in (a) for two different powers n. Near u = both functions are similar but for higher values they differ dramatically 
since the case n = 100 has an abrupt drop close to u — K. Their corresponding time evolution is plotted in (b) at a time at = 5 
starting from the same initial condition Ksech(x/a) where a is an arbitrary length. 



that would require a growth term bilinear rather than linear (see for example Ref. 22] ) . The generalization of the 
nonlinearity llOfl to a nonlinear growth term can be expressed as 

/(«) = {u/K) m [1 - {u/K) n ] (13) 

where m + 1 is the first power in u in a Taylor expansion near zero. Prescription J2J gives 

G(u) = — L ' , J . (14) 

where %Fi is the hypergeometric function. In general, Eq. (|14|l cannot be inverted and G _1 has to be calculated 
implicitly. However, analytical considerations regarding the shape of tails and shoulders of traveling fronts can be 
easily obtained through the help of Eq. l|14f) as will be seen below in Sec II I II 



E. Nonlinearity with provisional saturation 

The previous nonlinearities we have discussed above increase (decrease) monotonically from one zero to a maximum 
(minimum) and then monotonically decrease (increase) to the next zero. Our prescription J2J permits the analysis 
of more interesting situations that may arise if the mutual struggle for resources among individuals constituting the 
species leads to saturation of the population size that is not permanent but only provisional. Further increase in the 
population density in such a situation might lead to a returned increase in the nonlinearity, followed by an eventual 
decrease. Such a state of affairs is analogous to the case analyzed in subsection B above but with an important 
difference. In the present case, uf(u) does not dip below zero. As a consequence, there are no multiple zeros at u =/= 0. 
In Fig.Q we present, as an example of such a nonlinearity, a polynomial of the form, 

„ /w=9( „ W -i51(^ + M_n^ (16) 

and the corresponding time evolution. As time evolves, the curvature of the solution reflects the specific particulars 
of the "reaction" term. In the solution the regions with u values close to the top of the reaction function have 
more inclination. To study the effects of different dips, in Fig.JSJ three similar uf(u) are shown which differ in the 
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FIG. 4: Example of a nonlinearity with provisional saturation. In (a) the reaction function uf(u)/K of the polynomial 1151 is 
presented, while in (b) the corresponding evolution for an initial condition is shown. 



'depth' of their valleys. Their construction is straightforward by means of the Lagrange interpolating polynomial, 
i.e. by taking a polynomial of fourth order and requiring it to have specific values at the five points desired. In 
the example of the figure, four of the five points that are required to be common are: (0; 0), (1; 0), (1/4; 1), (3/4; 1) 
described in (u/K; uf(u)) plane. The other point will characterize the dip, and for the example shown in the figures 
are (1/2; 9/10), (1/2; 1/2), (1/2; 1/10). From the same initial condition, the evolution will differ as shown in (b): the 
more pronounced the dip in the nonlinearity, the flatter the solution close to the center of the dip. 




FIG. 5: Dependence of the depth of the dip for nonlinearities with provisional saturation. In (a) three different nonlinearities 
are plotted. The axis are interchanged respect to the previous figure. This has been done to express clearly the correspondence 
between the nonlinearities in (a) and the evolution in (b). 
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III. TRAVELING FRONT SOLUTIONS: ANALYTICAL EXPRESSIONS FOR SPATIAL LIMITS 



The analysis in the previous section has focused on the initial value problem. Much interesting information in 
nonlinear problems may be extracted, however, directly from traveling waves (see, e.g., extensive reports in the 
soliton literature such as in |2^|). The G- function appearing in our prescription Q can be used conveniently, and 
directly, to find traveling wave solutions. As will be seen below, there are situations in which this technique is useful 
even when explicit analytic expressions cannot be found for the initial value problem. Imposing the traveling wave 
ansatz u{x, t) = U{z — x — ct) in Eq.J2J), 

-(3^ = aUf(U), (16) 

where 



(3 = c-v (17) 

is the difference between the traveling front velocity and the medium velocity. Integrating i|16|) in z, leads to the 
solution 

G[U p (z)] = aL (18) 
P 

where G is defined in Eq. (0J. If the function G _1 can be written analytically, the front solution can be obtained in 
explicit form: 

(19) 

Otherwise it is defined implicitly through Eq. 1|18|) . 

Equation i|19|) shows that, for a given nonlinearity, the steepness of the front depends on the propagation velocity 
c while its shape is determined by the function G . In order to show how, for a given propagation velocity c, the 
shape of the front changes for various nonlinearities (according to the simple relation 119| l). we provide Fig. (jHJ). 



U p [z) = G~ 



P. 




FIG. 6: Example of several traveling waves of the same speed, (3 = 1, but for different nonlinearity functions /. The solid line 
is the generalized logistic function, ufi(u) = u(l — u 2 ), while the dotted line comes from the first half period of the sinusoidal 
function, 11/2(11) = sinljiu). 

We now show how one may find analytic expressions for the tail and the shoulder of the traveling wave through 
inspection of the limiting values of the nonlinearity auf(u), respectively, as u — > + and u — * K~ , where and K are 
both zeros of the nonlinearity. The basic idea is to apply to Eq. I|18|) a judiciously chosen function h, multiply by 
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U or (K — U) as appropriate, and then calculate the relevant limit of the product. Thus, if we are interested in the 
spatial dependence of the tail of the solution, we seek to obtain 

lim U(z)h T (G(U)) = lim U(z)h T (az/ [3) = C T . (20) 

C/->0+ z//3->+oo 

The function Ht is chosen in such a way that the limit for U — > + of the left hand side of the equation equals a 
non zero constant depicted above by Ct- If the shoulder spatial dependence is sought, we apply another function hs, 
and multiply by (K — U) rather than by U before taking the limit U — > K~, the choice of hs being such that the 
limit is a constant Cs- 

The spatial dependence of the traveling solution in the two extreme limits is found to be 

lim U(z) ~ - Cr (21) 
z//3^+oo h T (az/f3) 

for the tail and 

lim U(z) ~K - - f S ... (22) 

for the shoulder. 

The steepness of the front at the inflection point is another important feature of the traveling front that can be 
obtained from the shape of the nonlinearity even if the exact analytic form is not known. From Eq. H16(l we can 
simply relate 

dU(z) JP(z)] 

— = - aU{z) ^^- (23) 

which clearly shows the relation between the maximum of f(u) and the front derivative at the inflection point. 
Comparing such values with those cases of f(u) in which the traveling front shape are known analytically, allows to 
determine further qualitative characteristics of the front shapes. We apply these various general considerations to two 
specific cases below. 



A. Linear growth with power saturation 

From the general expression © , it is straightforward to write the analytic form of the traveling fronts 

m*) = ~ ^TAT- 



l)e an #V 



where we have chosen to center the traveling front in such a way that Up{z — 0) is half the saturation value of U. 
Application of the limiting procedure described above shows that regardless of the value of the power saturation n 
the front tail is given by 

lim U p {z) cx e - az/l3 . (25) 

The value n does not play a role since the tail shape is determined by the form of the nonlinearity close to U = 
where saturation effects are negligible. On the other hand, the value n does play a role in determining the shoulder 
shape. The limiting procedure gives in fact 

lim U (z) ocK- —e-^'P. (26) 

which coincides with the one calculated directly from Eq. Q24JI. 

It is interesting to notice that in the limiting case corresponding to n — » oo the traveling front solution can be 
written as 

U M = \Ke-'"> ifz>0 W 
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B. Nonlinear growth 

For the case in which the growth near u ~ is bilinear as in l|13|) for m = 2, the prescription J2J is given in terms 
of hypergeometric functions (see Ref. |24|). whose inverses are not known. However our limit prescriptions i|21[l and 
(|22[1 do provide analytical expressions for the spatial limits. For [/«0 the following behavior is obtained 

G{u) oc r ^ — — - . (28) 
v ' Km(U/K) m y ' 

The function Ht of Eq. (|20|l which makes the left limit a constant is found to be 

h T (G) = G^. (29) 

Therefore, in the tail, 

lim U(z)~ J— . (30) 

z/p^+oc {z/(3) — 

Indeed, we can determine the shape of the tail (shoulder) for any function / that has a Taylor expansion near [7 = 
(U = K). 



IV. INCORPORATION OF DIFFUSION ALONG WITH CONVECTION 



The analytic treatment we have given above of the reaction-convection problem represented by Eq. J2J) has been 
made possible by the fact that the partial differential equation considered is of first order. Once the transformation is 
made to convert (J5J exactly into its linear counterpart, the method of characteristics |25j, well-known in the context 
of linear equations, is what is behind the analysis we have presented. If diffusion is added to @, individual points 
of the u(x) curve at a given time no longer evolve independently as they do for systems characterized by (0 but 
influence one another in the evolution. Finding analytic solutions in the manner detailed here and in Refs. |l5l Il7| 
is then impossible. There is no doubt from the practical viewpoint that the case D — is important to study. This 
is so because negligible diffusion does occur in several physical situations as in the case of bacterial dynamics when 
the environment or genetic engineering can be made to result in convection being far more important than diffusion. 
Nevertheless it is important to ask whether the introduction of diffusion can be considered a perturbation or whether 
it alters the problem drastically, even when it is relatively small in magnitude. If the latter were true, the analysis 
presented in the previous sections would be uninteresting for all realistic cases in which finite diffusion exists. We 
therefore analyze the effects of adding diffusion to the nonlinear convective problem: 

du du „, . ft 2 u 

- +v - = auf(u) + D w . (31) 

We do this in three parts, always focusing attention on traveling wave solutions for simplicity. In the first part 
(subsection A below), we obtain exact analytical expressions for traveling wave solutions of the reaction diffusion 
problem represented by equation (|31|l when / is a piecewise linear function. Also, we compare the results between the 
absence of diffusion to those in its presence (arbitrary amount of diffusion). In the second part, when the amount of 
diffusion is small, we identify an appropriate small parameter and develop a perturbation scheme. In the third part, 
we apply the perturbation treatment to one of the interesting nonlinearities (sinusoidal) introduced at the beginning 
of this paper. Furthermore, we analyze its validity by comparing it to the exact solution of l|31|) when the nonlinearity 
is piecewise linear. The last two parts make up subsection B below. 



A. Diffusion effects treated through a piecewise linear representation of the nonlinearity 

Although exact analytic solutions of the nonlinear equation we consider are impossible in the presence of diffusion, 
they are possible for traveling waves if we employ a piecewise linearization of the nonlinearity. The piecewise linear 
approximation has been used earlier in the analysis of the Fisher equation for the study of the effects of transport 
memory [9j as well as pattern formation |26(. In both cases it was possible to extract useful information through 
this resource, whose advantage is that once the approximation of representing the given nonlinearity by the piecewise 
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linear form is made, no further approximation is invoked. Let us begin with a piecewise representation of the logistic 
nonlinearity through 



f(u)={ M/^ (32) 




in which the first piece of uf(u) connects the unstable zero at u = to the maximum at u = aK, while the second 
piece joins the maximum with the stable zero at u = K . The parameter < a < 1 in (|32|l represents the relative 
position of the maximum of the nonlinearity between the stable and unstable zeros of the system. If a — 1/2, the 
derivatives of the logistic nonlinearity and its piecewise representation become equal at u — and u = K . Since the 
tail of a front solution is determined by the limit l|21[l . the choice a — 1/2 ensures that the spatial dependence of the 
tail (shoulder) front in the logistic nonlinearity coincides with the tail (shoulder) front when f{u) is given by l|32|) . In 
the analytic treatment below, for the sake of generality, we keep the parameter a unspecified. 

Following the method of Ref. Q it is straightforward to determine the exact form of the traveling wave solution of 
Eq. (|31l) with / given by 132|) . Such a solution, denoted as U(z, 7), where 7 = Vf//3 and Vf = 2\/Da is given by 

U(az/j3 <0,-y)=K-K(l- a) exp > ( 33a ) 

U(az/}3 > 0, 7) = B 2 exp + B 3 exp (-^As) , (33b) 

with 



X 1 = (2/ 7 2 ) (J 1 + —72 -11, B 2 = (1/2)K (a-0, (33c) 

A 2 = (2/ 7 2 ) (l + v/Tf) , B s = (l/2)K(a + 0, (33d) 

. J I + j^j 2 , 

As = (2/ 7 2 ) (l - , £ = -(l-a) V ° +^=^- (33e) 

By requiring that it is always positive or zero, a condition for (|33|l emerges as 

H < !• (34) 

This bound represents a constraint between the Fisher speed w / and the relative (to the medium) speed of the front 
solution (3. Equation (|34|) is equivalent to the minimum speed requirement for traveling wave of reaction-diffusion 
systems and it expresses the fact that the front speed must be faster or equal to the diffusive speed. 

The parameter 7 is a measure of the relative diffusion in the system. In Fig. Q we show how the front solution 
changes as diffusion increases. As in reaction diffusion equations, it is clear in Fig. Q)) that, as diffusion increases, 
the tail of U(z,-f) acquires a shallower profile. The steepness of the tail as function of 7 can be determined from Eq. 
(|33b(l . The front tails are sums of two exponentials each one with characteristic lengths (A2 and A3) and multiplicative 
weighing factors (B 2 and B3), all functions of 7. When 7 = 0, B 2 = and \ 2 — 0. Eq. (|33b|l then contains only one 
exponential term. The tail of the front solution is thus given by 

U{az/f3 > 0,0) = B 3 exp =aK exp(-az//3). (35) 



/3A 3 

At the other extreme, i.e., when 7 = 1, the characteristic lengths X 2 and A3 are equal to each other, the divergent 
part of the coefficients, £, cancels in the sum, and the tail of the front solution becomes 

U (az/(3 > 0, 1) - (B 2 + B 3 )cxp(-^J+ = aK exp (~az/0\ . (36) 

For values < 7 < 1, both exponentials contribute to the tail shape but for small 7 values B 2 is negligible and 
A3 ^S> X 2 therefore making the second exponential in 133bl) dominate. This dependence is plotted in Fig. |SJ 

The use of the piecewise representation has allowed us to conduct an exact analysis above. When such a repre- 
sentation is not used, a perturbation scheme must be employed. The behavior shown in Fig. [S] lends support to the 
idea that a quantity directly related to 7 should serve as the small parameter on which to base the development of a 
perturbation technique for traveling front solutions. We develop such a scheme in the next subsection. 
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FIG. 7: Comparison of the shape of the traveling wave solution l|33[l with and without diffusion. On the left, in (a), a front 
solution with 7 = and one with 7 = 0.9 are depicted. On the right, in (b), AU(z, 7) = U(z,j) — 17(2,0) between Eq. (13311 
with 7 = 0.1, 0.5, 0.9 and the case with 7 = are shown. 




0.0 0.5 1.0 

7 

FIG. 8: Plot of the dependence on 7 of the parameters associated with the traveling front solution l|33b(l . The characteristic 
lengths of the exponentials A2 and A3 are depicted between the values 7 = and 7 = 1. In the inset the multiplicative factors 
B2 and Bz are also plotted. 



B. Perturbation treatment for traveling fronts 



The analysis presented here is a generalization of the technique proposed by Canosa 27] for studying traveling 
fronts of the Fisher equation. 
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In the presence of diffusion, the traveling fronts of Eq. I|31[) are solutions of the following equation 

,,2 



7 



U (( + U C + Uf(U) = 0, 



(37) 



where £ = azj (3 and z — x — ct where c is the absolute speed of the traveling front. 

The quantity £ = (j/2) 2 < 1/4 is the perturbation parameter we will choose. A standard expansion of the traveling 
front in powers of £ in Eq. l|31|l gives the identity 



n=l 



= 0, 



(38) 



where F (u) = uf(u). By Taylor-expanding the function F(u), and rearranging each term in powers of £, it is possible 
to find the approximation of order n in £ to the traveling front solution of Eq. I|31|l . The zeroth order of such an 
expansion is given by the solution of 



which has been analyzed in the previous sections. The next order is given by 

dUi cPUo p, rrn _ n 

-dt + ^e +UlF{Ua) -°-> 

whose solution is @ 

C/ 1 (C)=-^(C)ln(4|^(C)|) 



The next order, the second, is given by 



- } F"{U )Ul = 



which is a linear equation for U2, and all the other functions are given by solutions of the previous orders. 
To arbitrary order n, we have 



+ d ~^r + UnFW + h u u *<»» = 



(39) 

(40) 
(41) 
(42) 

(43) 



where h{\ Uj<r n } are the rest of the terms that come from the nth power of £. 

Canosa|23 and Murray^ have provided up to to first order for the Fisher case. For the trigonometric nonlinearity 
with f(u) given by (|1(J|) studied earlier in the present paper, application of this perturbation techniques gives for the 
first order correction 



Ui(z) — — sech(az/j3) In (sech(az / (3)) . 

7T 



(44) 



To further analyze the validity of the perturbation scheme, we use the exact expression for the traveling fronts 
of the piecewise linear function and compare it with the results obtained from the zeroth, first, and second order 
perturbation. Fig. [S] shows the differences between the solutions for two different values of the diffusion constant 
represented by 7 — 0.1 and 7 = 0.9. We see that, for small az//3, the lower orders are more accurate than the higher 
orders. In space this corresponds to that part of the front profile close to the value where U(z) = 1/2. However, for 
the rest of the azj 'ft range, the higher orders in the perturbation get closer to the actual solution. It is worth noticing 
that, for 7 = 0.9, which represents a high value of diffusion, the error committed is still low. Perhaps unexpectedly, 
for low values of az/0, the lower orders are more accurate. The convergence of the orders is fast: for 7 = 0.1 the 
second order is just overlaped by the first order. 

The above analysis has uncovered a perhaps obvious but important point regarding the perturbation parameter. 
In H31|) . there are two characteristic velocities: the Fisher velocity vj = \/Da and the medium velocity v. It may be 
tempting to suppose that the appropriate perturbation parameter is their ratio |15| . This is not correct because the 
entire effect of the medium velocity v can be eliminated by a transformation of the frame of reference. Our analysis 
above shows that it is crucial to take into consideration the third velocity that exists in the traveling wave context in 
this problem, viz. the velocity of the traveling wave c. The appropriate perturbation parameter when we consider a 
traveling wave is vt/(c — v) and our conclusion is that for diffusion to be considered as a perturbation on the pure 
convective analysis that we have presented in this paper, it is necessary that the velocity of the traveling wave exceed 
the medium velocity by an amount larger than the Fisher velocity. 
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FIG. 9: Difference between exact solution of the piecewise linear function and its zeroth, first and second order approximations 
in our perturbative treatment for two different values: 7 = 0.1 and 7 = 0.9 in (a) and (b) respectively. Notice that, although 
not shown in the graph, as the order of the perturbation increases the curve difference converges rapidly to a curve which 
represents the error upper bound. 



V. SUMMARY 



Several physically realizable situations exist in various areas of nonlinear population dynamics wherein the behavior 
is controlled primarily by the nonlinearities and the convective element of the motion, the diffusive contribution to 
the motion being small. The typical example is the behavior of bacteria in a petri dish [H Q E|, where the 
so-called "wind" motion induced by moving masks is made stronger than the inherent diffusive motion. A study of 
such situations has been provided in the present paper by first analyzing the case for no diffusion (sections IIIIIII|) 
and then extending the analysis to small diffusion (section ITvTl . Explicit results have been obtained for a variety of 
nonlinearities on the basis of a straightforward prescription [13|, for the full initial value problem. Traveling waves 
have been analyzed explicitly for several cases. Interesting consequences such as the "pyramid effect" for sinusoidal 
nonlinearities have been shown to occur. The effect of small diffusion has been incorporated through exact analysis 
for piecewise representation of the nonlinearities on the one hand and through perturbative calculations on the other. 
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